function [yyyy , qtr] = gen_yq_vec(start_y , start_q , l_q)
    if start_q > 4
        error('Starting quarter cannot be larger than 4')
    end
    
    % Start with start year and first quarter
    l_q_full = l_q + start_q - 1;
    yyyy_full = [kron((start_y : start_y + floor(l_q_full/4) - 1)' , ones(4 , 1));...
                (start_y + floor(l_q_full/4)) * ones(mod(l_q_full , 4) , 1)];
            
    qtr_full =  [kron(ones(floor(l_q_full/4) , 1) , (1 : 4)');...
                (1 : mod(l_q_full , 4))'];
            
    % Make adjustment for starting quarter        
    yyyy = yyyy_full(start_q : end);  
    qtr = qtr_full(start_q : end);  
end